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We study the statistics of ecosystems with a variable number of co-evolving species. The species 
interact in two ways: by prey-predator relationships and by direct competition with similar kinds. 
The interaction coefficients change slowly through successful adaptations and speciations. We treat 
them as quenched random variables. These interactions determine long-term topological features of 
the species network, which are found to agree with those of biological systems. 
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Population dynamics is a classical subject of evolution- 
ary biology. The mutually dependent dynamics of two or 
more populations or species is often described by coupled 
differential equations governing the relative change of the 
population sizes Ni{t), 
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The interaction coefficients gij can represent a prey- 
predator relationship [gij < 0, gji > 0), direct compe- 
tition (gij < 0, gji < 0), or mutualism [gij > 0, gji > 0) 
between species i and j, and the terms hi denote intrinsic 
production or death rates. These so-called Lotka- Volterra 
equations, as well as many generalizations thereof, have 
been used to model coexistence, invasions, and adaptive 
change of populations. Of great importance is their con- 
ceptual connection to mathematical game theory A 
set of populations A^i , . . . , Ns represents a mixed strat- 
egy. For given interactions gij , an optimal strategy - 
called Nash equilibrium - can often be realized as a stable 
fixed point , . . . , Ng of an associated Lotka- Volterra 
dynamics. This explains how strategic optimization is 
reached in biological systems through reproductive suc- 
cess, with no need for rational thinking. 

These equilibria determine the species' fate. For a 
given set of equations (Q), a species is viable if A^* > 
and becomes extinct if A^* = 0. Even the viable species 
are not perennial, however. Successful adaptations, mi- 
grations, and speciations (the splittings of a single species 
into a pair) eventually change the number of players as 
well as the rules of the game, i.e., the couplings gij. On 
large time scales, this dynamics can be quite intermittent. 
Correlated extinctions and speciations alternate with pe- 
riods of relative stasis, leading to large fluctuations in the 
number of species |^ . Little is known on how this long- 
term behavior is connected to the underlying interactions 
between species in (m). 
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FIG. 1. (a) The Pamlico estuary foodweb in North Car- 
olina, consisting of 14 species (filled circles) at four trophic 
levels. Detritus, dinolagellates and diatoms are at the bot- 
tom level [l = 1) and feed from external resources (empty 
symbols). There is a single trophic group at the highest level 
(i = 4), formed by the predatory fishes Roccus and Cynoscion. 
Arrows point from prey to predator; dashed lines connect 
species pairs with a nonzero link overlap (see text). Data 
from Cohen (1990). (b) Average species numbers for a set 
of natural ecosystems, taken from Cohen (1990) (empty sym- 
bols) and Rosenzweig (1995) (filled circles). This last case 
corresponds to an average over 61 independent food webs, 
most of which are empty at high levels. 

The best studied natural ecosystems are food webs, i.e., 
communities of animal species in a closed environment 
where food chains can be observed. Fig. 1(a) shows the 
graph of such a network, each arrow representing a prey- 
predator relationship. Despite large variations in size and 
environmental conditions, large ecosystems share a few 
important topological characteristics: (i) Every species 
lives at a certain trophic level, which can be defined as the 
minimum length of its relevant 'downward' food chains. 
Species at level one feed from external resources, (ii) The 
number of trophic levels is small, typically between three 
and seven, (iii) Most species have a small number of rel- 
evant prey species (typically around three) , mainly from 
the next lower level, (iv) The number of species at level I 
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increases with / for lower values of / and decreases again 
sharply for higher I ^j^, see Fig. 1(b). Networks of co- 
evolving species thus have a characteristic shape. 

This remarkable structure calls for a theoretical expla- 
nation. The classical work on Lotka-Volterra equations 
has established stability criteria for networks with ran- 
dom interactions gij Q . In a real ecosystem, however, the 
interactions are not random, but are themselves subject 
to selection. Recently, Lotka-Volterra systems coupled to 
speciation and immigration dynamics have been 
studied by numerical simulations, and food web struc- 
tures have indeed been found. Another class of models 
focuses directly on the dynamics of extinctions and spe- 
ciations. These models have no explicit population dy- 
namics and mostly random topology, with the important 
exception of Ref. . 

In this Letter, we present elements of a statistical the- 
ory for large ecosystems, using concepts and methods 
of theoretical physics. We discuss the population dy- 
namics with the minimal species interactions consistent 
with the observed complexity of ecosystems. These inter- 
actions are prey-predator relationships^ which establish 
a flux of biomass between species, and direct competi- 
tion between similar species, which leads to their mutual 
exclusion from ecological niches. The interaction coeffi- 
cients gij in (|l|) are modeled as random variables that 
change through successful mutations. We focus on the 
(often realistic) case that these mutations are sufficiently 
rare so that the populations can reach stable equilibria 
in between. In the language of statistical physics, the 
species interactions are quenched random variables on 
the time scales of population dynamics. The statistics 
of such ecosystems is thus governed by a quenched dis- 
tribution of Nash equilibria. This distribution in turn 
emerges from a long-term balance between adaptations, 
speciations, and extinctions. The topology of the result- 
ing networks is found to be closely related to the un- 
derlying dynamics of co-evolution. In the following, we 
concentrate on generic topological features amenable to 
an approximate analytical treatment; in particular, we 
derive the shape of ecosystems. A detailed analysis of 
structure and dynamics of these networks will be pub- 
lished elsewhere |l0| ]. 

To describe generic features of ecological networks, we 
choose the simplest population dynamics containing pre- 
dation and direct competition. The interaction matrix 
in (|^) is decomposed accordingly, gij = 'jij — f3ij. Pre- 
dation is parameterized by the constants 7,,, — 7+ if 
j is prey of i and 7^^ — —7- if i is prey of j, with 
< 7+ < 7-. Competition takes place for nesting places, 
mating opportunities, and other resources not explicitly 
represented in the model. It is strongest between individ- 
uals of the same species, but also occurs between different 
species that interfere in each other's livelihood We 
set f3ii = 1 (this normalization amounts to an appropri- 
ate choice of the time scale in (ft)) and /3y = (3pij for 



i ^ j, with < P < 1. The link overlap pij measures 
the degree of competition between the species. It is de- 
fined as Pij = Cij/ ^CiCj, where Ci, Cj are the number of 
predatory links of i, j, and Cy is the number of com- 
mon predatory links. (Species pairs with a nonzero link 
overlap are connected by dashed lines in Fig. 1.) Fur- 
thermore, all species are assigned a uniform death rate 
hi = —a. The external resources are represented as a 
small number of extra 'populations' Ni with hi — 7+i? 
and predators only (i.e., 7ij < and /3y = for all j). 

With these interactions, the fixed point populations 
given by (|l]) can be written in the form 

N* = P,-Q, + h,, (2) 

where P, = 7+ Ejex(») - 7- Ejen(.) N* is the pro- 
ductivity of species i from predation (with irii) the set 
of its prey and n(i) the set of its predators) and Qi — 
PJ2j=ii Pij-^j is its competition load. Furthermore, we 
require a minimum population size Nc <^ R for viable 
species, and count all species with N* < Nc as extinct. 
Indeed, natural populations are known to be unstable 
under short-term environmental fluctuations or adverse 
mutations if they are too small or too dilute |l^ . 

Of course, an ecosystem is not determined by its popu- 
lation dynamics alone but also by the long-term processes 
of successful mutations, in particular, speciations [f7|,p^. 
In this model, a mutation is represented as a stochastic 
change of predation links that is consistent with exist- 
ing food chains. It turns out that details of this pro- 
cess are not relevant for our present purpose of deriv- 
ing global network characteristics. It is sufficient to as- 
sume that speciations and adaptations maintain a broad 
distribution of productivities Pi, and hence, of popula- 
tion sizes N* (in a sense made precise below). This is 
well supported by field observations and by our numer- 
ics jl^,|ljl . An increase in the number of species reduces 
the average productivity and increases the average com- 
petition load. Hence, such an ecosystem admits only a 
certain number of viable species, whose productivities 
satisfy Pi > Qi -\- a -\- Nc- The number of these ecological 
niches depends on the interaction parameters /3, 7+, 7_, 
and on the dimensionless ratios R/Nc, a/Nc- Once the 
niches are filled, ongoing speciations and the subsequent 
adaptations reshuffle the productivities and the popu- 
lation numbers N* of all the species, forcing the least 
viable ones into extinction. On large time scales, this 
is a stationary stochastic process. The relative success 
of an individual species keeps changing as a result of its 
own adaptations and those of the other species, resulting 
in a constant threat of extinction called the Red Queen 
effect [|l^ . The shape of these mature networks is deter- 
mined essentially by the distribution of ecological niches. 
To see this, consider first two cases of simple networks 
with fixed topology. 

1. A single food chain is a community of L species on 
L trophic levels. The species at level one feeds from an 



2 



external resource, the species at level I from that at level 
I — 1 {I = 2, 3, . . . , L). The productivities of this chain 
are given by the equations 



(/ = 1,...,L) 



(3) 



and Pq — —jiN^, with the boundary condition Nl_^^ ~ 
0. They determine directly the population numbers 
TV;* = Pi — a since all competition loads vanish. The 
entire chain is viable if Pi > Pc for all species /, with the 
minimum productivity 
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The equations (||) can be solved exactly by recursion 
starting from the top level I — L. For the biologically 
important case of small 7+, we find that the maximum 
value of L compatible with (||) is 



L = 



log7^ 



log 



R 



const. a + Nc 



- 1 + 0(74 
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by applying the condition (^) at the top level. 

The parameters a and Nc are seen to be equivalent 
viability cutoffs for the chain since they reduce primarily 
the top population iV£. More generally, the population 
numbers are found to be rapidly decreasing with in- 
creasing I for all relevant parameter values. Hence, as 
observed in nature, viable chains are always short p6[ |. 

2. A single trophic level is a group of S species that 
may have a significant overlap in their predation links 
and a resulting competition load. First we consider the 
productivities Pi as fixed by the interactions with other 
trophic levels and concentrate on the effects of the direct 
competition terms Qi. In a 'mean field' approximation, 
we replace the individual link overlaps by an expectation 
value p depending on the predation clusters. In the sim- 
plest case of random predation, Eq. (^) then determines 
the fixed point populations 
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the average N = S ^J2i=i given by 
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The viability of all species {N* > Nc) again sets a mini- 
mum productivity 



Pc^a+{1- Pp)Nc + PpSN 



(8) 



We now use the assumption that the productivities Pi 
are drawn from a broad probability distribution given 
by <^{q) = Prob(Pi/P < q). (The qualitative results 
do not depend strongly on the form of $((7); here we 
use a simple approximation p7[|.) The species commu- 
nity becomes unstable if the least viable species has a 



productivity below Pc- The number of species in a ma- 
ture trophic level can therefore be estimated from the 
relation S'^{Pc/P) = 0(1). Eq. (||) then becomes an im- 
plicit relation for S" as a function of P/Nc, a/Nc, and the 
average pairwise competition load (3p. Consider, for ex- 
ample, a trophic level with random predation from a set 
of S' prey species from the levels below. Using a simple 
approximation for the average link overlap p{S, S') [ p^ , 
it can be shown that the solution of Eq. (||) always sat- 
isfies S < max(a(/3)S"/c, 1), where a{(3) > 1 and c is 
the average number of prey species per predator species. 
That is, competition determines the number of ecologi- 
cal niches in a trophic level as a function of the prey di- 
versity and the competition strength /3. For sufficiently 
large /?, only non-overlapping species can coexist, i.e., 
S — max (S"/c, 1). This result generalizes the well known 
theorem of competitive exclusion ||l9|| , which states the 
condition for coexistence of two competing species. Note 
that this limiting effect on the number of species exists 
independently of the population numbers. It is indeed 
crucial for the buildup of high population numbers at the 
lower trophic levels. For example, a trophic level feed- 
ing from effective resources of size R ^ a, Nc acquires 
an extensive population number per species N R/S, 
while S is asymptotically independent of R. Without 
competitive exclusion (/3 = 0), speciations would further 
increase 5". This leads eventually to an extensive number 
of marginally viable species, i.e., S ~ R/N with N of or- 
der Nc- Such a level could not support sizeable predation 
from above. 

We now turn to a full ecological network with L trophic 
levels. In the mean field approximation, we treat all 
species at the same level on an equal footing and derive 
self-consistent equations for the level averages of popula- 
tion and species number, Ni and Si {I = 1, . . . , L). The 
average productivities Pi satisfy the recursion relations 



Pi = 7+c^;-i - -/-ciSi+i/Si)Ni+i 



(9) 



where we assume that the species at every level pre- 
date randomly on the species at the next lower level. 
The average number c of predation links per predator is 
taken to be independent of Z; this is indeed suggested by 
field data. The average number of predators per prey 
is then simply cSi+i/Si. The productivity Pi is linked 
to Ni and Si as in (0), using for p{Si,Si-i) the same 
approximation as above ||l^. Hence, the relations (^) de- 
termine the population numbers given the species num- 
bers. The latter are again limited by the stability crite- 
ria Si^{Pc^i/ Pi) = 0(1) with the minimum productivi- 
ties Pc,; given as in (^); these relations determine the Si 
given the Ni . The coupled set of equations can be solved 
iteratively. Finally, the number of levels L follows from 
the condition N^ ~ iVcj which is equivalent to S'l « 1. 

Over a wide range of relevant parameters, these net- 
works have the characteristic shape shown in the example 
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of Fig. 2: The species numbers Si increase with I at low 
levels due to the increasing prey diversity, which opens 
up more and more niches. They reach a maximum at 
an intermediate level and decrease again at higher levels, 
because more and more species have population numbers 
too low to support further predation. Hence, these two 
regimes reflect the two kinds of species interactions. The 
population numbers show an approximately exponential 
decrease in both regimes, just like for a single vertical 
chain. Hence, L is always small, in agreement with ob- 
servations and with the results of The functional 
form of the patterns Si , Ni can be described by analytical 
approximations depending on the parameter values. 




FIG. 2. The shape of ecosystems, (a) The species numbers 
Si {1 < I < L) for networks with L trophic levels. The param- 
eters are c = 3, 7+ = 0.3, 7- = 2.0, P — 0.2, qo = 0.35 (see 
[17]), a/Nc = 1, and R/Nc = 2x 10^ 10*, 4 x lO"*, 2 x 10^ 
and 5 x 10^ for the cases L = 4, 5, 6, 7, and 8, respectively, 
(b) The average population numbers Ni for the same cases as 
in (a). 

Species networks are thus quite far from randomly con- 
nected. Their topological shape is dynamically generated 
by the coupled evolution of populations and the slower 
adaptative changes. The ubiquity of this shape suggests 
that predation and competition of similar species are the 
fundamental interactions governing the long-term coevo- 
lution of large ecosystems. They are remarkably simple. 
Predation is the basic transport of energy in the system, 
competition forces the species into states with little over- 
lap. In physics, mutual avoidance is a well known prop- 
erty of fermions. Competitive exclusion may thus be re- 
garded as the Pauli principle of co-evolution: It generates 
the complexity of species networks just as its quantum- 
mechanical counterpart does for atoms and molecules. 

We have discussed here the global shape of these net- 
works. It remains a challenging task to explore the con- 
nection between dynamics and topology locally, that is, 
at the level of individual species and their genealogies. 

M.L. is grateful to the MPI for Colloids and Interfaces 
for the kind hospitality throughout this work. 
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